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Abstract 

The cutting plane approach to optimal matchings has been discussed by several authors over 
the past decades [HI HH [TTl [25l [10] , and its convergence has been an open question. We prove 
that the cutting plane approach using Edmonds' blossom inequalities converges in polynomial 
time for the minimum-cost perfect matching problem. Our main insight is an LP-based method 
to retain/drop candidate cutting planes. This cut retention procedure leads to a sequence of 
intermediate linear programs with a linear number of constraints whose optima are half-integral 
and supported by a disjoint union of odd cycles and edges. This structural property of the 
optima is instrumental in finding violated blossom inequalities (cuts) in linear time. Further, 
the number of cycles in the support of the half-integral optima acts as a potential function to 
show efficient convergence to an integral solution. 

1 Introduction 

Integer programming is a powerful and widely used approach to modeling and solving discrete 
optimization problems [19, 22]. Not surprisingly, it is NP-complete and the fastest known algorithms 
are exponential in the number of variables (roughly [16J). In spite of this intractability, 

integer programs of considerable sizes are routinely solved in practice. A popular approach is the 
cutting plane method, proposed by Dantzig, Fulkerson and Johnson [7] and pioneered by Gomory 
|11 ] 112 ] [T3]. This approach can be summarized as follows: 

1. Solve a linear programming relaxation (LP) of the given integer program (IP) to obtain a 
basic optimal solution x. 

2. If x is integral, terminate. If x is not integral, find a linear inequality that is valid for the 
convex hull of all integer solutions but violated by x. 

3. Add the inequality to the current LP, possibly drop some other inequalities and solve the 
resulting LP to obtain a basic optimal solution x. Go back to Step 2. 

"This work was supported in part by NSF award AF0915903. This work was done while the second author was 
affiliated with the College of Computing, Georgia Institute of Technology, supported by NSF Grant CCF-0914732. 
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For the method to be efficient, we require the following: (a) an efficient procedure for finding a 
violated inequality (called a cutting plane) , (b) convergence of the method to an integral solution 
using the efficient cut-generation procedure and (c) a bound on the number of iterations to conver- 
gence. Gomory gave the first efficient cut-generation procedure and showed that the the cutting 
plane method implemented using his cut-generation procedure converges to an integral solution [13J. 
There is a rich theory on the choice of cutting planes, both in general and for specific problems 
of interest. This theory includes interesting families of cutting planes with efficient cut-generation 
procedures [ITJ (H HI El G3 [2Q1 13 EEl [23], valid inequalities, closure properties and a classification 
of the strength of inequalities based on their rank with respect to cut-generating procedures [6] 
(e.g., the Chvatal-Gomory rank [2]), and testifies to the power and generality of the cutting plane 
method. 

To our knowledge, however, there are no polynomial bounds on the number of iterations to 
convergence of the cutting plane method even for specific problems using specific cut-generation 
procedures. The best bound for general 0-1 integer programs remains Gomory 's bound of 2 n |13j . 
It is possible that such a bound can be significantly improved for IPs with small Chvatal-Gomory 
rank [10J. A more realistic possibility is that the approach is provably efficient for combinatorial 
optimization problems that are known to be solvable in polynomial time. An ideal candidate could 
be a problem that (a) has a polynomial-size IP-description (the LP-relaxation is polynomial-size), 
and (b) the convex- hull of integer solutions has a polynomial-time separation oracle. Such a problem 
admits a polynomial-time algorithm via the Ellipsoid method [15] . Perhaps the first such interesting 
problem is minimum-cost perfect matching: given a graph with costs on the edges, find a perfect 
matching of minimum total cost. 

A polyhedral characterization of the matching problem was discovered by Edmonds [8] : Basic 
solutions of the following linear program (extreme points of the polytope) correspond to perfect 
matchings of the graph. 



The relaxation with only the degree and nonnegativity constraints, known as the bipartite relaxation, 
suffices to characterize the convex-hull of perfect matchings in bipartite graphs, and serves as a 
natural starting relaxation. The inequalities corresponding to sets of odd cardinality greater than 1 
are called blossom inequalities. These inequalities have Chvatal rank 1, i.e., applying one round of 
all possible Gomory cuts to the bipartite relaxation suffices to recover the perfect matching polytope 
of any graph [4J. Moreover, although the number of blossom inequalities is exponential in the size 
of the graph, for any point not in the perfect matching polytope, a violated (blossom) inequality 
can be found in polynomial time [21j . This suggests a natural cutting plane algorithm (Figure [T]), 
proposed by Padberg and Rao [21] and discussed by Lovasz and Plummer in their classic book on 
matching theory [17J. Experimental evidence suggesting that this method converges quickly was 
given by Grotschel and Holland |14] , by Trick [25J , and by Fischetti and Lodi |10| . It has been open 
to rigorously explain their findings. In this paper, we address the question of whether the method 
can be implemented to converge in polynomial time. 
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The known polynomial-time algorithms for minimum-cost perfect matching are variants of Ed- 
monds' weighted matching algorithm [8]. It is perhaps tempting to interpret the latter as a cutting 
plane algorithm, by adding cuts corresponding to the shrunk sets in the iterations of Edmonds' 
algorithm. However, there is no correspondence between the solution x of the LP given by non- 
negativity and degree constraints and a family T of blossom inequalities, and the partial matching 
M in the iteration of Edmonds' algorithm when T is the set of shrunk nodes. In particular, the 
next odd set S shrunk by Edmonds' algorithm might not even be a cut for x (i.e., x(5(S)) > 1). 
It is even possible, that the bipartite relaxation already has an integer optimal solution, whereas 
Edmonds' algorithm proceeds by shrinking and unshrinking a long sequence of odd sets. 



1. Start with the bipartite relaxation. 

2. While the current solution is fractional, 

(a) Find a violated blossom inequality and add it to the LP. 

(b) Solve the new LP. 



Figure 1: Cutting plane method for matchings 

The bipartite relaxation has the nice property that any basic solution is half-integral and its 
support is a disjoint union of edges and odd cycles. This makes it particularly easy to find violated 
blossom inequalities - any odd component of the support gives one. This is also the simplest 
heuristic that is employed in the implementations |14| [25] for finding violated blossom inequalities. 
However, if we have a fractional solution in a later phase, there is no guarantee that we can find 
an odd connected component whose blossom inequality is violated, and therefore sophisticated and 
significantly slower separation methods are needed for finding cutting planes, e.g., the Padberg-Rao 
procedure [21]. Thus, it is natural to wonder if there is a choice of cutting planes that maintains 
half-integrality of intermediate LP optimal solutions. 




Graph G The starting optimum xq Basic feasible solution obtained 

and the cut to be imposed after imposing the cut 

Figure 2: Counterexample to the half-integrality conjecture. All edge costs are one. 

At first sight, maintaining half- integrality seems to be impossible. Figure [2] shows an example 
where the starting solution consists of two odd cycles. There is only one reasonable way to impose 
cuts, and it leads to a non half-integral basic feasible solution. Observe however, that in the 
example, the bipartite relaxation also has an integer optimal solution. The problem here is the 
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existence of multiple basic optimal solutions. To avoid such degeneracy, we will ensure that all 
linear systems that we encounter have unique optimal solutions. 

This uniqueness is achieved by the simple deterministic perturbation of the integer cost function, 
which increases the input size polynomially. We observe that this perturbation is only a first step 
towards maintaining half-integrality of intermediate LP optima. More careful cut retention and 
cut addition procedures are needed to maintain half-integrality. 

1.1 Main result 

To state our main result, we first recall the definition of a laminar family: A family T of subsets 
of V is called laminar, if for any X, Y G J 7 , one of X n Y = 0, X C Y, Y C X holds. 

Next we define a perturbation to the cost function that will help avoid some degeneracies. Given 
an integer cost function c : E — > Z on the edges of a graph G = (V,E), let us define the perturbation 
c by ordering the edges arbitrarily, and increasing the cost of edge i by 1/2*. 

We are now ready to state our main theorem. 

Theorem 1.1. Let G = (V, E) be a graph on n nodes with edge costs c : E — > Z and let c denote 
the perturbation of c. Then, there exists an implementation of the cutting plane method that finds 
the minimum c-cost perfect matching such that 

(i) every intermediate LP is defined by the bipartite relaxation constraints and a collection of 
blossom inequalities corresponding to a laminar family of odd subsets, 

(ii) every intermediate LP optimum is unique, half-integral, and supported by a disjoint union of 
edges and odd cycles and 

(Hi) the total number of iterations to arrive at a minimum c-cost perfect matching is 0(n logn). 

Moreover, the collection of blossom inequalities used at each step can be identified by solving an LP 
of the same size as the current LP. The minimum c-cost perfect matching is also a minimum c-cost 
perfect matching. 

To our knowledge, this is the first polynomial bound on the convergence of a cutting plane 
method for matchings. It is easy to verify that for an n- vertex graph, a laminar family of nontrivial 
odd sets may have at most n/2 members, hence every intermediate LP has at most 3n/2 inequalities 
apart from the non-negativity constraints. 

1.2 Cut selection via dual values 

The main difficulty that we have to overcome is the following: ensuring unique optimal solutions 
does not suffice to maintain half-integrality of optimal solutions upon adding any sequence of blos- 
som inequalities. In fact, even a laminar family of blossom inequalities is insufficient to guarantee 
the nice structural property on the intermediate LP optimum. Thus a careful choice is to be made 
while choosing new cuts and it is also crucial that we eliminate certain older ones that are no longer 
useful. 

At any iteration, inequalities that are tight for the current optimal solution are natural candi- 
dates for retaining in the next iteration while the new inequalities are determined by odd cycles in 
the support of the current optimal solution. However, it turns out that keeping all tight inequal- 
ities does not maintain half-integrality. Our main algorithmic insight is that the choice of cuts 
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for the next iteration can be determined by examining optimal dual solutions to the current LP 
- we retain those cuts whose dual values are strictly positive. Since there could be multiple dual 
optimal solutions, we use a restricted type of dual optimal solution (later called positively-critical 
dual in this paper) that can be computed either by solving a single LP of the same complexity or 
combinatorially. Moreover, we also ensure that the set of cuts imposed in any LP are laminar and 
correspond to blossom inequalities. 

Eliminating cutting planes that have zero dual values in any later iteration is common in most 
implementations of the cutting plane algorithm; although this is done mainly to keep the number 
of inequalities from blowing up, another justification is that a cut with zero dual value is not a 
facet contributing to the current LP optimum. 

1.3 Algorithm C-P-Matching 

Let G = (V, E) be a graph, c : E — > R a cost function on the edges, and assume G has a perfect 
matching. The primal bipartite relaxation polytope and its dual are specified as follows. 

min c{uv)x{uv) (Pq(G,c)) 

uv£E 

x(5(u)) = l \/u GV max^Tr(u) (D (G,c)) 

x > 

7T{U) + TT[y) < c[uv) VUV G E 



We call a vector x G M. E proper-half-integral if x(e) G {0, 1/2, 1} for every e G E and supp(x) is 
a disjoint union of edges and odd cycles. The bipartite relaxation of any graph has the following 
well-known property. 



Proposition 1.2. Every basic feasible solution x of Pq(G,c) is proper-half-integral. □ 

Let O be the set of all odd subsets of V of size at least 3, and let V denote the set of one element 
subsets of V . For a family of odd sets JCO, consider the following pair of linear programs. 

min c{uv)x{uv) (Pj^(G,c)) 

uv£E 

x(5{u)) = l VuGV max S n(S) ( D AG,c)) 



x(s(S)) > i ys g t 

y n(5) < c(uv) Vuv G E 

x > ^— ' 



E 

SeVuF:uve5(s) 



n(5) > V5 G J 7 



For J 7 = 0, Pjr(G,c) is identical to Po(G, c) whereas for J 7 = 0, it is identical to |PJ). Every 



intermediate LP in our cutting plane algorithm will be \Pj^(G, c)| for some laminar family P. We 
will use n(u) to denote n({u}) for dual solutions. 
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Assume we are given a dual feasible solution T to Djr(G, c). We say that a dual optimal solution 
II to Djr(G,c) is T- extremal, if it minimizes 



h(U,T)= 



\U(S)-T(S)\ 



among all dual optimal solutions II. A T-extremal dual optimal solution can be found by solving 
a single LP if we are provided with the primal optimal solution to Pj^{G, cj| (see Section 5.3). 



The cutting plane implementation that we propose is shown in Figure [3j From the previous 
set of cuts, we retain only those which have a positive value in an extremal dual optimal solution; 
let %' denote this set of cuts. The new set of cuts %" correspond to odd cycles in the support of 
the current solution. However, in order to maintain laminarity of the cut family, we do not add 
the vertex sets of these cycles but instead their union with all the sets in H 1 that they intersect. 
We will show that these unions are also odd sets and thus give blossom inequalities. In the first 
iteration, there is no need to solve the dual LP as T will be empty. 



1. Let c denote the cost function on edges after perturbation (i.e., after ordering the 
edges arbitrarily and increasing the cost of edge i by 1/2*). 

2. Starting LP. Let T be the empty set, so that the starting LP, Pjr(G,c), is the 
bipartite relaxation and the starting dual V is identically zero. 

3. Repeat until x is integral: 

(a) Solve LP. Find an optimal solution x to Pj^(G, c). 

(b) Choose old cutting planes. Find a T-extremal dual optimal solution II to 
D T (G,c). Let 

%> = {S g T : n(5) > 0}. 

(c) Choose new cutting planes. Let C denote the set of odd cycles in supp(x). 
For each C G C, define C as the union of V{C) and the maximal sets of H' 
intersecting it. Let 

n" = {C:C € C}. 

(d) Set the next T = H! U H" and r = IL 

4. Return the minimum-cost perfect matching x. 



Figure 3: Algorithm C-P-Matching 



1.4 Overview of the analysis 

Our analysis to show half-integrality is based on the intimate relationship that matchings have with 
factor-critical graphs (deleting any node leaves the graph with a perfect matching): for example, 
contracted sets are factor-critical both in the unweighted [9| and weighted |8j matching algorithms 
by Edmonds. We define a notion of factor-criticality for weighted graphs that also takes a laminar 



6 



odd family T into account. This notion will be crucial to guarantee half-integrality of intermediate 
solutions. 

We use the number of odd cycles in the support of an optimal half-integral solution as a potential 
function to show convergence. We first show odd(xj+i) < odd(xj), where Xi,Xi+\ are consecutive 
LP optimal solutions, and odd(.) is the number of odd cycles in the support. We further show that 
the cuts added in iterations where odd(xj) does not decrease continue to be retained until odd(xj) 
decreases. Since the maximum size of a laminar family of nontrivial odd sets is n/2, we get a bound 
of 0(n log n) on the number of iterations. 

The proof of the potential function behavior is quite intricate. It proceeds by designing a half- 
integral version of Edmonds primal-dual algorithm for minimum-cost perfect matching, and arguing 
that the optimal solution to the extremal dual LP must correspond to the one found by this primal- 
dual algorithm. We emphasize that this algorithm is used only in the analysis. Nevertheless, it 
is rather remarkable that even for analyzing the cutting plane approach, comparison with a new 
extension of Edmonds' classic algorithm provides the answer. 



2 Factor-critical sets 



In what follows, we formulate a notion of factor-critical sets and factor-critical duals, that play a 
central role in the analysis of our algorithm and are extensions of concepts central to the analysis 
of Edmonds' algorithm. 

Let H = (V, E) be a graph and J 7 be a laminar family of subsets of V . We say that an edge set 
M C E is an T '-matching, if it is a matching, and for any S £ J-, \M D S(S)\ < 1. For a set S QV, 
we call a set M of edges to be an (S, J 7 ) -perfect-matching if it is an J-"- matching covering precisely 
the vertex set S. 

A set S G J- is defined to be (H , J-) -factor- critical or T -factor- critical in H, if for every node 
u G S, there exists an (S \ {u}, J r )-perfect-matching using the edges of H. For a laminar family 
T and a feasible solution II to Djr(G,c), let Gn = (V,Eji) denote the graph of tight edges. For 
simplicity we will say that a set S G T is (II, ^-factor-critical if it is (Gn, .F)-factor critical, 
i.e., S is J 7 - factor-critical in Gri- For a vertex u G S, corresponding matching M u is called the 
LT- critical-matching for u. If T is clear from the context, then we simply say S is II-factor-critical. 

A feasible solution LT to Djr(G, c) is an J- '-critical dual, if every S G J- is (II, J 7 )-factor-critical, 
and II(T) > for every non-maximal set T of T. A family T C O is called a critical family, if T is 
laminar, and there exists an ^-critical dual solution. This will be a crucial notion: the set of cuts 
imposed in every iteration of the cutting plane algorithm will be a critical family. The following 
observation provides some context and motivation for these definitions. 

Proposition 2.1. Let T be the set of contracted sets at some stage of Edmonds' matching algo- 
rithm. Then the corresponding dual solution IT in the algorithm is an J- '-critical dual. □ 



We call IT to be an T -positively-critical dual, if II is a feasible solution to \Djr(G, c) and every 



S G T such that 11(5) > is (LT, ^-factor-critical. Clearly, every ^-critical dual is also an 
J-"-positively-critical dual, but the converse is not true. The extremal dual optimal solutions found 
in every iteration of Algorithm C-P-Matching will be J-"-positively-critical, where T is the family 
of blossom inequalities imposed in that iteration. 

The next lemma summarizes elementary properties of LT-critical matchings. 
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Lemma 2.2. Let J 7 be a laminar odd family, H be a feasible solution to Djr(G, c) , and S £ J- be a 



(II, J 7 ) -factor- critical set. For u^v G S, let M u , M v be the II- critical-matching s for u, v respectively, 
(i) For every T G J 7 such that T C S, 



\M u nS(T)\ 



ifueS\T, 
ifu G T. 



(ii) Assume the symmetric difference of M u and M v contains an even cycle C. Then the sym- 
metric difference M U AC is also a H-critical matching for u. 

Proof, (i) M u is a perfect matching of S \ {u}, hence for every T C S, 

\M u r\5(T)\ = \T\{u}\ (mod 2). 

By definition of M u , \M U n S(T)\ < 1 for any T C S, T G J 7 , implying the claim. 

(ii) Let M' = M U AC. First observe that since C is an even cycle, u, v G" V{C). Hence M' is a 
perfect matching on using only tight edges w.r.t. II. It remains to show that \M' n<5(T)| < 1 

for every T G J 7 , T C S. Let 7„ and j v denote the number of edges in C n 5(T) belonging to M u 
and M v , respectively. Since these are critical matchings, we have ^y u ,7v < 1- On the other hand, 
since C is a cycle, |C n 5(T)\ is even and hence 7„ + j v = \C D 5(T)\ is even. These imply that 
7„ = 7„. The claim follows since \M' n <5(T")| = |M U n <5(T)| - 7„ + 7«. □ 

The following uniqueness property is used to guarantee the existence of a proper-half-integral 
solution in each step. We require that the cost function c : E — > K satisfies: 



For every critical family J 7 , Pj^(G, c) has a unique optimal solution. (*) 



The next lemma shows that an arbitrary integer cost function can be perturbed to satisfy this 
property. The proof of the lemma is presented in Section [7j 

Lemma 2.3. Let c : E — ^ Z be an integer cost function, and c be its perturbation. Then c satisfies 
the uniqueness property (pp. 



3 Analysis outline and proof of the main theorem 



The proof of our main theorem is established in two parts. In the first part, we show that half- 
integrality of the intermediate primal optimum solutions is guaranteed by the existence of an 



J-"-positively-critical dual optimal solution to Djr(G, c 



Lemma 3.1. Let J 7 be a laminar odd family and assume Pj^(G,c) has a unique optimal solution 



x. If there exists an J 7 -positively- critical dual optimal solution, then x is proper-half -integral. 



Lemma 3.1 is shown using a basic contraction operation. Given II, an ^-positively-critical 



dual optimal solution for the laminar odd family J 7 , contracting every set S G J 7 with H(S) > 
preserves primal and dual optimal solutions (similar to Edmonds' primal-dual algorithm). This 

Moreover, if we had a unique primal optimal solution x to Pj^(G,c 



is shown in Lemma 4.1 



its 



image x' in the contracted graph is the unique optimal solution; if x' is proper-half-integral, then 
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so is x. Lemma 3.1 follows: we contract all maximal sets S G J- with H(S) > 0. The image x' 
of the unique optimal solution x is the unique optimal solution to the bipartite relaxation in the 
contracted graph, and consequently, half-integral. 

Such ^-positively-critical dual optimal solutions are hence quite helpful, but their existence 
is far from obvious. We next show that if J 7 is a critical family, then the extremal dual optimal 
solutions found in the algorithm are in fact ^-positively-critical dual optimal solutions. 

Lemma 3.2. Suppose that in an iteration of Algorithm C-P-Matching, T is a critical family with T 
being an T '-critical dual solution. Then a V -extremal dual optimal solution II is an T ' -positively- critical 
dual optimal solution. Moreover, the next set of cuts W = %' U %" is a critical family with II being 
an T~L- critical dual. 

To show that a critical family T always admits an ^-positively-critical dual optimum, and 
that every extremal dual solution satisfies this property, we need a deeper understanding of the 
structure of dual optimal solutions. Section [5] is dedicated to this analysis. Let T be an ^-critical 



dual solution, and II be an arbitrary dual optimal solution to Dj?(G,c). Lemma 5.1 shows the 
following relation between II and T inside sets S G J- that are tight for a primal optimal solution 
x: Let Ts(u) and Us(u) denote the sum of the dual values of sets containing u that are strictly 
contained inside S in solutions Y and II respectively, and let A = max ng s(rs(it) — lis (it)). Then, 
every edge in supp(x) n 6(S) is incident to some node u G S such that Tg(u) — lis (it) = A. Also 



(Lemma 5.8), if S G J- is both T- and IT-factor-critical, then T and II are identical inside S. 

If 11(5) > but S is not LT-factor-critical, the above property (called consistency later) en- 
ables us to modify II by moving towards V inside S, and decreasing 11(5') so that optimality is 
maintained. Thus, we either get that LT and T are identical inside S thereby making S to be 
LT-factor-critical or H(S) = 0. A sequence of such operations converts an arbitrary dual optimal 
solution to an ^-positively-critical dual optimal one, leading to a combinatorial procedure to ob- 



tain positively-critical dual optimal solutions (Section 5.2). Moreover, such operations decrease the 
secondary objective value h(IL, T) and thus show that every T-extremal dual optimum is also an 
^-positively-critical dual optimum. 



Lemmas |3.1| and 3.2 together guarantee that the unique primal optimal solutions obtained 
during the execution of the algorithm are proper-half-integral. In the second part of the proof of 
Theorem |1.1| we show convergence by considering the number of odd cycles, odd(x), in the support 
of the current primal optimal solution x. 

Lemma 3.3. Assume the cost function c satisfies (Qj. Then odd(x) is non-increasing during the 
execution of Algorithm C-P-Matching. 



We observe that similar to Lemma |3.1[ the above Lemma 3.3 is also true if we choose an 
arbitrary ^-positively-critical dual optimal solution II in each iteration of the algorithm. To show 
that the number of cycles cannot remain the same and has to strictly decrease within a polynomial 
number of iterations, we need the more specific choice of extremal duals. 

Lemma 3.4. Assume the cost function c satisfies Q) and that odd(x) does not decrease between 
iterations i and j, for some i < j. Let Tk be the set of blossom inequalities imposed in the k'th 
iteration and H!' k = Fk \ Fk-i be the subset of new inequalities in this iteration. Then, 

U 

k=i+l 
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We prove this progress by coupling intermediate primal and dual solutions with the solutions 
of a Half-integral Matching algorithm, a variation of Edmonds' primal-dual weighted matching 
algorithm that we design for this purpose. 

We briefly describe this next, assuming familiarity with Edmonds' algorithm [8]. Our argument 



needs one phase of this algorithm and this is what we analyze in detail in Section 6.1 Although 



we introduce it for analysis, we note that the algorithm can be extended to a strongly-polynomial 
combinatorial algorithm for minimum-cost perfect matching. Unlike Edmonds' algorithm, which 
maintains an integral matching and extends the matching to cover all vertices, our algorithm 
maintains a proper half-integral solution. 

The half-integral matching algorithm starts from a partial matching x in G, leaving a set W of 
nodes exposed, and a dual II whose support is a laminar family VU J- with T C 0; x and II satisfy 
primal-dual slackness conditions. The algorithm transforms x to a proper-half-integral perfect 
matching and II to a dual solution with support contained in VU/C, satisfying complementary 



slackness. Before outlining the algorithm, we show how we apply it to prove Lemmas 3.3 and 3.4 



Let us consider two consecutive primal solutions Xi and x«+i in the cutting plane algorithm, with 
duals Hi and IL + i. We contract every set S £ O with IIj + i(S') > 0; let G be the resulting graph. 



By Lemma 3.1 the image x' i+1 of Xj+i is the unique optimal solution to the bipartite relaxation 
in G. The image x\ of Xi is proper-half-integral in G with some exposed nodes W; let 11^ be the 
image of LTj. Every exposed node in W corresponds to a cycle in supp(xj). We start in G with 
the solutions x\ and ITj, and we prove that it must terminate with the primal solution x' i+1 . The 
analysis of the half-integral matching algorithm reveals that the total number of exposed nodes and 
odd cycles does not increase; this will imply Lemma |3.3| 

To prove Lemma 3.4, we show that if the number of cycles does not decrease between phases i 
and i + 1, then the algorithm also terminates with the extremal dual optimal solution B^ , x . This 
enables us to couple the performance of Half-integral Matching between phases i and i + 1 and 
between i + 1 and i + 2: the (alternating forest) structure built in the former iteration carries over 
to the latter one. As a consequence, all cuts added in iteration i will be imposed in all subsequent 
phases until the number of odd cycles decreases. 

Let us now turn to the description of the Half-integral Matching algorithm. In every step of the 
algorithm, we maintain z to be a proper-half-integral partial matching with exposed set T C W, A 
to be a dual solution satisfying complementary slackness with z, and the support of A is V U C for 
some C C T . We start with z = x,T = W,A = U and C = T '. We work on the graph G* resulting 
from the contraction of the maximal sets of C 

By changing the dual A, we grow an alternating forest of tight edges in G* , using only edges e 
with z{e) = or z(e) = 1. The forest is rooted at the exposed set of nodes T. The solution z will 
be changed according to the following three scenarios, (a) If we find an alternating path between 
two exposed nodes, we change z by alternating along this path as in Edmonds' algorithm, (b) If 
we find an alternating path P from an exposed node to a 1/2-cycle C in supp(z), we change z by 
alternating along P, and replacing it by a blossom (an alternating odd cycle) on C. (c) If we find 
an alternating path P from an exposed node to a blossom C, then we change z by alternating along 
P and replacing the blossom by a 1/2-cycle on C. If none of these cases apply, we change the dual 
A in order to extend the alternating forest. If A(S) decreases to for some S S C, then we remove 
S from C and unshrink it in G* . 

The modifications are illustrated on Figure |4j Note that in case (c), Edmonds' algorithm would 
instead contract C. In contrast, we do not perform any contractions, but allow 1/2-cycles in the 
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half-integral matching algorithm. For the special starting solution x = 0, II = 0, our half-integral 
matching algorithm will return a proper-half-integral optimum to the bipartite relaxation. 



o 



o 



1 



Figure 4: The possible modifications in the Half-integral Matching algorithm. 
The main theorem can now be proved using the above lemmas. 

Proof of Theorem \1.1\ We use Algorithm C-P-Matching given in Figure [3] for a perturbed cost 
function. By Lemma |2.3[ this satisfies Q. Let i denote the index of the iteration. We prove by 
induction on i that every intermediate solution x% is proper-half-integral and (i) follows immediately 
by the choice of the algorithm. The proper-half-integral property holds for the initial solution xq by 



Proposition 1.2. The induction step follows by Lemmas 3.1 and 3.2 and the uniqueness property. 



Further, by Lemma 3.3, the number of odd cycles in the support does not increase. 



Assume the number of cycles in the i'th phase is £, and we have the same number of odd cycles 
£ in a later iteration j. Between iterations i and j, the set T~L'^ always contains I cuts, and thus 
the number of cuts added is at least £(j — i). By Lemma 3.4, all cuts in (Jfe=i+i ^1 are imposed 
in the family J~j+i- Since J~j+i is a laminar odd family, it can contain at most n/2 subsets, and 
therefore j — i < n/2£. Consequently, the number of cycles must decrease from £ to £ — 1 within 
n/2£ iterations. Since odd(a;o) < n/3, the number of iterations is at most O(nlogn). 

Finally, we show that optimal solution returned by the algorithm using c is also optimal for 
the original cost function. Let M be the optimal matching returned by c, and assume for a 
contradiction that there exists a different perfect matching M' with c(M') < c(M). Since c is 
integral, it means c(M') < c(M) — 1. In the perturbation, since c(e) < c(e) for every e G E, we 
have c(M) < c(M) and since ]C eg£ ;(c(e) - c(e)) < 1, we have c(M') < c(M') + 1. This gives 
c(M') < c(M') + 1 < c(M) < c(M), a contradiction to the optimality of M for c. □ 



4 Contractions and half-integrality 

We define an important contraction operation and derive some fundamental properties. Let J 7 be a 
laminar odd family, LT be a feasible solution to Djr(G,c), and let S £ J- be a (LI, ^-factor-critical 
set. Let us define 

n s («) := Yl n ( T ) 

TeVUF:TCS,«eT 

to be the total dual contribution of sets inside S containing u. 

By contracting S w.r.t. LI, we mean the following: Let G' = (V',E') be the contracted graph 
on node set V = (V \ S) U {s}, s representing the contraction of S. Let V denote the set of 
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one-element subsets of V . For a set T C V, let T" denote its contracted image. Let T' be the set 
of nonsingular images of the sets of J-, that is, T' £ T' if T 6 J 7 , and T' \ {s} 7^ 0. Let 22' contain 
all edges uv £ E with u,v ^ S and for every edge with u £ S, v £ V — S add an edge st> . Let 
us define the image II' of II to be n'(T') = II(T) for every T' £ V U J 7 ' and the image x' of x to be 
x'iu'v 1 ) = x(uv). Define the new edge costs 




c(uv) iiuv£E[V\S], 
c(uv) - U s {u) if u £ S, v £ V \ S. 




Figure 5: Contraction operation: image of dual and new cost function 



Lemma 4.1. Let J 7 be a laminar odd family, x be an optimal solution to Pjr(G,c), IT be a feasible 
solution to Djr(G,c). Let S £ J- be a (II, F) -factor- critical set, and let G',c',J-' denote the graph, 
costs and laminar family respectively obtained by contracting S w.r.t. U and let x',IL' be the images 
of x,II respectively. Then the following hold. 

(i) IT' is a feasible solution to Dj?i{G' ,c'). Furthermore, if a set T£j r ,T\S^(/)is (II, J-)- 
f actor- critical, then its image T is (II' , T') -factor- critical. 



(ii) Suppose II is an optimal solution to Dj^(G, c) and x(S(S)) 
to Pj?'(G',c') and IT is optimal to Djr>(G' ,c'). 



1. Then x' is an optimal solution 



(Hi) If x is the unique optimum to Pjr(G,c), and Ii is an optimal solution to Djr(G,c), then x' 
is the unique optimum to Pjn{G' ,c'). Moreover, x' is proper-half-integral if and only if x is 
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proper-half-integral. Further, assume C" is an odd cycle in supp(x') and let T be the pre-image 
of V(C) in G. Then, supp(x) inside T consists of an odd cycle and matching edges. 

Proof, (i) For feasibility, it is sufficient to verify 

n'(T') < c'(u'v') VuV g e'. 

T'eV'LIT':u'v'e5(T') 

If u, v 7^ s, this is immediate from feasibility of II to Dp. (G, c). Consider an edge sv' G E{G'). Let 
uv be the pre-image of this edge. 



£ n'(r') = n(5) + Yl n ( T ) <cM-n 5 («) 

T'eV'UF'-.sv'eSiT') TeJ":™e<5(T),T\5^0 



c'(sv') 



We also observe that u'v' is tight in G' w.r.t II' if and only if the pre-image uv is tight in G w.r.t 

n. 

Let T G J- be (II, J r )-factor-critical with T \ 5 7^ 0. It is sufficient to verify that T is (II', J-')- 
factor-critical whenever T contains 5. Let v! G T 1 be the image of u G V, and consider the 
projection M' of the II-critical-matching M u . The support of this projection is contained in the 
tight edges with respect to II'. Let Z' C T' , Z' G J 7 ' and let Z be the pre-image of Z' . If u' ^ s, 



then \M' n 5(Z')| = |M U n 5(Z)| < 1 and since \M U n 5(5) | = 1 by Lemma 2.2, the matching 



M' is a (T' \ {u'}, J"')-perfect-matching. If u' = s, then u £ S. By Lemma [2^2} M u n 5(5) = 
and hence, M' misses s. Also, \M U n 5(Z)| < 1 implies |M' n 5(Z')\ < 1 and hence M' is a 
(T' \ {s}, ^^-perfect-matching. 

(ii) Since x(5(S)) = 1, we have x'(S(v)) = 1 for every v G V. It is straightforward to verify 
that x'{5(T')) > 1 for every T G T' with equality if x(5(T)) = 1. Thus, is feasible to P^{G f , d). 
Optimality follows as x' and II' satisfy complementary slackness, using that the image of tight 
edges is tight, as shown by the argument for part (i). 

(Hi) For uniqueness, consider an arbitrary optimal solution y' to Pj:t(G',d). Let M u be the 
LT-critical matching for u in 5. Define a u = y'(5(u)) for every u G 5, i.e., a u is the total y' value 
on the images of edges uv G E with v G V — 5. Note that X]ue5 a « = v'i^i^)) = 1- Take 

jV(uV) if £ E \ E[S] , 
I w(uv) if G £7[5]. 

Then, y is a feasible solution to Pj^(G,c) and y satisfies complementary slackness with II. Hence, 
y is an optimum to Pp(G, c) and thus by uniqueness, we get y = x. Consequently, y' = x' . 

The above argument also shows that x must be identical to w inside 5. Suppose x' is proper- 
half-integral. First, assume s is covered by a matching edge in x' . Then a u = 1 for some u G 5 
and a v = for every v ^ u. Consequently, u; = M u is a perfect matching on 5 — n. Next, 
assume s is incident to an odd cycle in x' . Then a Ul = a U2 = 1/2 for some nodes ui,u% G 5, and 
w = h(M Ul + M U2 ). The uniqueness of x implies the uniqueness of both M Ul and M U2 . Then by 



y(uv) 



Lemma |2.2[ ii), the symmetric difference of M U1 and M U2 may not contain any even cycles. Hence, 
supp(w) contains an even path between u\ and U2, and some matching edges. Consequently, x is 
proper-half-integral. The above argument immediately shows the following. 

Claim 4.2. Let C' be an odd (even) cycle such that x'(e) = 1/2 for every e G C' in supp(x') and 
let T be the pre-image of the set V(C") in G. Then, supp(x) n E[T] consists of an odd (even) cycle 
C and a (possibly empty) set M of edges such that x{e) = 1/2 V e G C and x(e) = 1 V e G M. □ 
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Next, we prove that if x is proper-half-integral, then so is x' . It is clear that x' being the image 
of x is half-integral. If x 1 is not proper-half-integral, then supp(x') contains an even 1/2-cycle, and 



thus by Claim 4.2 supp(x) must also contain an even cycle, contradicting that it was proper-half- 
integral. Finally, if C is an odd cycle in supp(x'), then Claim 4.2 provides the required structure 



for x inside T. We also need to show that its support does not contain an even cycle C . If not, 
then the same argument as above shows that the pre-image of C is an even cycle C in the support 
of x, a contradiction. □ 



Corollary 4.3. Assume x is the optimal solution to Pjr(G, c) and there exists an J- -positively- critical 
dual optimum II. Let G,c be the graph, and cost obtained by contracting all maximal sets S £ T 
with H(S) > w.r.t. U, and let x be the image of x in G. 

(i) x and II are the optimal solutions to the bipartite relaxation Po(G, c) and Dq(G, c) respectively. 

(ii) If x is the unique optimum to Pj?(G,c), then x is the unique optimum to Pq(G,c). If x, is 
proper-half-integral, then x is also proper-half-integral. 



Proof of Lemma \3. 1\ Let II be an ^-positively-critical dual optimum, and let x be the unique 
optimal solution to Pj^(G,c). Contract all maximal sets S £ J- with II (3) > 0, obtaining the 
graph G and cost c. Let x be the image of x in G. By Corollary 
to Pq(G, c). Bv Proposition 1.2 



By Proposition 
proper- half-integral. 



4.3 



x is proper-half-integral and hence by Corollary 



x is uniq ue optimum 
x is also 
□ 



4.3 



5 Structure of dual solutions 

In this section, we show two properties about positively-critical dual optimal solutions - (1) an 
optimum ^ to Dj^(G, c) can be transformed into an ^-positively-critical dual optimum (Section 



5.2) if J 7 is a critical family and (2) a T-extremal dual optimal solution to Djr(G,c) as obtained in 

we 



the algorithm is also an ^-positively-critical dual optimal solution (Section |5.3| ). In Section 5.1 
first show some lemmas characterizing arbitrary dual optimal solutions. 

5.1 Consistency of dual solutions 

Assume J- C Q is a critical family, with II being an .F-critical dual solution, and let ^ be an 



arbitrary dual optimal solution to Djr(G,c) Note that optimality of II is not assumed. Let x be 
an optimal solution to \Pjr(G, c)\ we do not assume uniqueness in this section. We shall describe 
structural properties of compared to IT; in particular, we show that if we contract a II-factor- 
critical set S, the images of x and ^ will be primal and dual optimal solutions in the contracted 
graph. 

Consider a set S € F. We say that the dual solutions II and VI/ are identical inside S, if 
n(T) = Vl/(T) for every set T C S, T £ T U V. We defined Hs(u) in the previous section; we also 
use this notation for ty, namely, let *&s( u ) '■= Y^teVut-tcs uet ^(T) for u £ S. Let us now define 

A n ,*(S) := max (U s (u) - V s {u)) . 

We say that ^ is consistent with IT inside S, if lis (it) ~^s( u ) = An,#(£) holds for every u £ S that 
is incident to an edge uv £ 6(S) D supp(x). The main goal of this section is to prove the following 
lemma. 
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Lemma 5.1. Let T CO be a critical family, with II being an T -critical dual solution and let Vl/ be 
an optimal solution to Dj^(G, cjj Let x be an optimal solution to \Pj^(G, c)\ Then ^> is consistent 



with II inside every set S £ J- such that x(5(S)) = 1. Further, Aji^(S) > for all such sets. 

Consistency is important as it enables us to preserve optimality when contracting a set S £ T 
w.r.t. II. Assume ^ is consistent with II inside S, and x(5(S)) = 1. Let us contract S w.r.t. II to 
obtain G' and c' as defined in Section [H Define 

if T' £ (J'UV')\W, 
-A n ,*(S) if T' = {s} 



*'(T') 




Lemma 5.2. Let T C O be a critical family, with II being an T '-critical dual solution and let \P be 
an optimal solution to Djr(G, c) . Let x be an optimal solution to Pp(G, c)j Suppose Vl/ is consistent 



with II inside S £ J- and x(5(S)) = 1. Let G' \d ,J-' denote the graph, costs and laminar family 
obtained by contraction. Then the image x' of x is an optimum to Pjri{G',c'), and \T/' (as defined 
above) is an optimum to Djn{G' ,d) . 

Proof. Feasibility of x' follows as in the proof of Lemma |4.l[ ii). For the feasibility of ^' , we have 
to verify ^T'eVuF'-mie^Cr') ^'(T') < c'(uv) for every edge uv £ E{G'). This follows immediately 
for every edge uv such that u,v ^ s since ^ is a feasible solution for Djr(G,c). Consider an edge 
uv £ E(G), u £ S. Let sv £ E(G') be the image of uv in G' , and let A = Ari^S 1 ). 

c(uv) > *(T) 

TeVUT:uv£8(T) 

TeT:uveS(T),T-S^9 

= M> S ( U ) + A + *'( T ')- 

T'eV'UF':sv£8{T') 

In the last equality, we used the definition ^'(s) = ^(S) — A. Therefore, using Hs(u) < ^(tz) + A, 
we obtain 

^'( T ') < c (uv) - V s (u) - A < c{uv) - U s (u) = c'(uv). (1) 

T'eV'UT':sv&&{T') 

Thus, is a feasible solution to Djn(G' , d). To show optimality, we verify complementary slackness 
for x' and ^' . If x'(uv) > for u,v ^ s, then x(uv) > 0. Thus, the tightness of the constraint 
for uv w.r.t. ^f' in Djn{G' ,d) follows from the tightness of the constraint w.r.t. \& in Djr(G,c). 
Suppose x'(sv) > for an edge sv £ E(G'). Let uv £ E(G) be the pre-image of sv for some u £ S. 
Then the tightness of the constraint follows since both the inequalities in ([I]) are tight - the first 
inequality is tight since uv is tight w.r.t. ^, and the second is tight since n,s(u) — ^s{u) = A(S) 
by the consistency property. Finally, if ^'(T 1 ) > for some T 1 £ J-', then \1>(T) > and hence 
x(S(T)) = 1 and therefore, x'{8(T')) = 1. □ 

Lemma 5.3. Let T be a critical family with II being an J- '-critical dual, and let x be an optimal 



solution to Pjr(G,c) If x(5(S)) = 1 for some S £ T , then x(5(T)) = 1 for every T C S, T £ T , 



and all edges in supp(x) n E[S] are tight w.r.t. II. 

Proof. Let a u = x(5(u, V \ S)) for each u £ S, and for each T C S, T £ J 7 , let a{T) = J2ueT Q " = 
x(S(T, V\S)). Note that a(S) = x(5(S)) = 1. Let us consider the following pair of linear programs. 
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min c{uv)z{uv) 

uv£E[T] 

z(5(u)) = 1 — a u Vu G 5 
z(S{T)) > 1 - a(T) VF C 5, T G J" 
z(uw) >OVraG F[S] 



(Ff[S]) 



max ^ (l-Q(T))r(T) (F^S]) 
^ T(T) < c{uv) \fuv G F[S*] 



Tcs,TeVuF 

uv£S(T) 



T(Z) > VZ C T, Z G F 



For a feasible solution z to Fj-JS'] , let x z denote the solution obtained by replacing x(uv) by 



z(uv) for edges uv inside 5, that is, 



x z (e) 




if e G 5(5')U£[F\5] ) 
if e G E[S}. 



Claim 5.4. T/ie restriction of x inside S is feasible to Pjr and for every feasible solution z to 
Pjr £> , x z is a feasible solution to Pj^(G, c) , Consequently, z is an optimal solution to Pp[S] if and 



only if x z is an optimal solution to Pjr(G,c) 



Proof. The first part is obvious. For feasibility of x z , if u ^ S then x z (u) = x(u) = 1. If u G S, 
then x z {u) = z(u) + x(5(u,V \ S)) = 1 — a u + a u = 1. Similarly, if T G F, T \ 5 ^ 0, then 
x 2 (<5(T)) = a?(T) > 1. If T C S, then x 2 (c5(T)) = z(£(F)) + x{5(T, V\S))>1 — a(T) + a(T) = 1. 
Optimality follows since c T x z = XmoeMS] c(-uv)z(ut;) + YluveE\E[S] c(uv)x(uv). □ 



Claim 5.5. Let LT denote the restriction of TI inside S, that is, n(T) 
T C. S. Then II is an optimal solution to Djr[S] 



II(T) /or every T G VUF, 



Proof. Since n is an F-critical dual and S G F, we have a II-critical-matching M u inside S for each 



u £ S. Let u) = X^eS The claim follows by showing that w is feasible to Pj?[S] and that w 

and n satisfy complementary slackness. 

The degree constraint w(5(u)) 
T G F, then w(5(T)) 



shows feasibility of II to Dj- [S] 



eS\T 



Oti, 



1 — a u is straightforward and using Lemma 2.2 if T C S 1 , 
= 1 — a(F). The feasibility of II to Djr(G,c) immediately 



Complementary slackness also follows since by definition, all M u 's use only tight edges w.r.t. II 
(equivalently, w.r.t. ft). Also, for every odd set T C S, T G F, we have that w(5(T)) = 1 — a(T) 
as verified above. Thus, all odd set constraints are tight in the primal. □ 



By Claim 5.4, the solution obtained by restricting x to E[S] must be optimal to Fj-JS 1 ] and 



thus satisfies complementary slackness with ft. Consequently, every edge in E[S] n supp(x) must 
be tight w.r.t. ft and equivalently w.r.t. II. By the definition of the F-critical-property, every 
non-maximal set T of F satisfies 11(F) > 0, and in particular, every set T C S, T G F satisfies 
fl(F) = 11(F) > 0. Thus, complementary slackness gives x(5(T)) = 1. 

□ 



We need one more claim to prove Lemma 5.1 
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Claim 5.6. Let S G J- be an inclusionwise minimal set of T . Let A and T be feasible solutions 
to 



Djr(G,c), and suppose S is (A, J 7 ) -factor- critical. Then A\p(S) = max ug s | Ag(it) — Ts(u) 



Further, if A\^(S) > 0, define 



A+ := {u G T : T(u) = A(u) + A A ,r(5)}, 
A- :={ueT: T(u) = A(u) - A A ,r(5)}. 



Then \A~\ > \A+\. 



Proof. Let A = max u6 s \As(u) — Ts(u)\, and define the sets A' and A + with A instead of AatOS')- 
Since S is (A, ^-factor-critical, for every a G S, there exists an (S \ {a}, J 7 ) perfect matching M a 
using only tight edges w.r.t. A, i.e., M a C {uv : A(u) + A(v) = c(uv)} by the minimality of S. 
Further, by feasibility of T, we have T(u) + T(v) < c{uv) on every uv G M a . Thus, if u G A + , 
then v G A~ for every uv G M a . Since A > 0, we have A + U A~ ^ and therefore A~ ^ 0, 
and consequently, A = A^r^S). Now pick a G A~ and consider M a . This perfect matching M a 
matches each node in A + to a node in A~ . Thus, |j4 _ | > \ A + \. □ 



Proof of Lemma 5.1, We prove this by induction on |V|, and subject to this, on \S\. Let S be a 



(II, J r )-factor-critical set. First, consider the case when S is an inclusion-wise minimal set. Then, 



U s (u) = U(u), = *(it) for every u G S. By Claim[K6j we have that A := A n ,*(5') > 0. We 

are done if A = 0. Otherwise, define the sets A~ and A + as in the claim using An,* (S). 

Now consider an edge uv G i?[S']nsupp(x). By complementary slackness, we have ^(u)+^>(v) = 
c(uv). By dual feasibility, we have H(u) + n(i>) < c(uv). Hence, if u G A~ , then v G A + . 
Consequently, we have 

\A~\ = x{5(u)) = x(5(A-,V\S)) + x(5(A-,A + )) 

u£A~ 

< 1 + ^2 x(5(u)) = 1 + < \A~\. 

u£A+ 

Thus, we should have equality throughout. Hence, x(5(A~ , V \ S)) = 1. This precisely means that 
^ is consistent with n inside S. 

Next, let S be a non-minimal set. Let T G J- be a maximal set strictly contained in S. By 



Lemma 5.3, x(5(T)) = 1, therefore the inductional claim holds for T: is consistent with n inside 



T, and A(T) = A n ,*(T) > 0. 

We contract T w.r.t. n and use Lemma 5.2 Let the image of the solutions x, H, and \& be x' , n' 
and respectively and the resulting graph be G' with cost function c' . Then x' and ty' are optimum 
to Pji{G' , c') and Djn(G' , c') respectively, and by Lemma |4.l[ i), n' is an ^'-critical dual. Let t be 
the image of T by the contraction. Now, consider the image S' of S in G' . Since G' is a smaller 
graph, it satisfies the induction hypothesis. Let A' = An',*' (5") in G' . By induction hypothesis, 
A' > 0. The following claim verifies consistency inside S and thus completes the proof. □ 

Claim 5.7. For every u G S, Us(u) — ^s{ u ) < Tl' s ,(u') — ^' s ,(u'), and equality holds if there exists 
an edge uv G S(S) n supp(x). Consequently, A' = A. 

Proof. Let u' denote the image of u. If v! ^ t, then U' s , (it') = lis (it) j ^'s> ( u ') = ^s(u) and therefore, 
n 5 (it) - = W s ,(u') - %,(u'). Assume u' = t, that is, u G T. Then H s (u) = U T (u) + n(T), 
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^s( u ) = ^t(u) + *&(T) by the maximal choice of T, and therefore 



U s (u) - = n T (u) - * r («) + n(T) - (T) 

< A(T) + n(T) - (T) 

= n'(t) - (Since W{t) = U{T), t>'(t) = *(T) - A(T)) 

= l4(t)-¥' s ,(t). (2) 

Assume now that there exists a uv £ 5(S) n supp(x). If u G T, then using the consistency inside T, 
we get U T (u) - tf T (tt) = A(T), and therefore (g gives U s (u) - = II' s ,(t) - \^,(t) = A'. □ 



Claim [5T6] can also be used to derive the following important property. 

Lemma 5.8. Given a laminar odd family T C O, let A and T be two dual feasible solutions to 
Djr(G,c). If a subset S £ T is both (A, J 7 ) -factor- critical and (T ', F)- factor- critical, then A and V 
are identical inside S. 

Proof. Consider a graph G = (V,E) with \V\ minimal, where the claim does not hold for some 
set S. Also, choose S to be the smallest counterexample in this graph. First, assume S G J- is a 



minimal set. Then consider Claim 5.6 for A and T and also by changing their roles, for V and A. 
If A and T are not identical inside S, then A = max ng ,5 \ As(u) — Ts(u)\ > 0. The sets A~ and A + 
for A and T become A + and A~ for T and A. Then \A~\ > > \A~\, a contradiction. 

Suppose now S contains T £ T . It is straightforward to see that T is also (A, ^-factor-critical 
and (r, ^-factor-critical by definition. Thus, by the minimal choice of the counterexample S, we 
have that A and T are identical inside T. Now, contract the set T w.r.t. A, or equivalently, w.r.t. 
r. Let A', r' denote the contracted solutions in G' , and let T' be the contraction of T . Then, by 
Lemma |4.l[ i), these two solutions are feasible to Dj?>(G' ,cf), and S' is both A'-factor-critical and 
r'-factor-critical. Now, A' and V are not identical inside S' , contradicting the minimal choice of G 
and S. □ 

5.2 Finding a positively-critical dual optimal solution 

Let J- C O be a critical family with II being an ^-critical dual. Let ^ be a dual optimum solution 
to Djr(G, c). Our goal is to satisfy the property that for every S & F, if ^f(S) > 0, then ^ and II 
are identical inside S. By Lemma |5.8[ it is equivalent to showing that is ^-positively-critical. 
We modify by the algorithm shown in Figure [6} 

The correctness of the algorithm follows by showing that the modified solution ^ is also dual 
optimal, and it is closer to IT. 

Lemma 5.9. Let T C O be a critical family with II being an J- -critical dual and let ^ be a dual 
optimum solution to Dj?(G,c). Suppose we consider a maximal set S such that U and are not 
identical inside S, take A = min{l, ^f(S)/ Au,^(S)} if An,* (-5) > and A = 1 if An,* (SO = and 
set ^S> as in Q). Then, ^ is also a dual optimal solution to Djr(G,c), and either 'f(S) = or II 
and ^ are identical inside S. 



Proof. Let x be an optimal solution to Pjr(G,c) Since \£(S) > 0, we have x(S(S)) = 1 and by 
Lemma 5.1, we have A = An,* (5) > 0. Now, the second conclusion is immediate from definition: 
if A = 1, then we have that II and ^ are identical inside S; if A < 1, then we have ^(S) = 0. For 
optimality, we show feasibility and verify the primal-dual slackness conditions. 
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1. Repeat while ^ is not ^-positively-critical dual. 

(a) Choose a maximal set S G J- with ^(S) > 0, such that II and \& are not identical inside 
S. 

(b) Set A := A n ,*(S). 

(c) Let A := min{l, 9(S) /A} if A > and A := 1 if A = 0. 

(d) Replace by the following \t. 



(1 - A)*(T) + An(T) if T C 5, 
t(T) := <j tf(S) - AA if T = S, 
Vl'(T) otherwise . 



(3) 



2. Return ^-positively-critical dual optimum ^. 



Figure 6: Algorithm Positively-critical-dual-opt 



The solution ^ might have positive components on some sets T C S,T G J 7 where n(T) > 0. 
However, x(5(T)) = 1 for all such sets by Lemma 5.3 The choice of A also guarantees ^f(S) > 0. 



1 for all such sets by Lemma ! 
We need to verify that all inequalities in Djr(G,c) are maintained and that all tight constraints in 
Djr(G, c) w.r.t. \I> are maintained. This trivially holds if uv G E[V \S]. If uv G E[S] \ supp(x), the 
corresponding inequality is satisfied by both LT and and hence also by their linear combinations. 
If uv G E[S] nsupp(x), then uv is tight for ^ by the optimality of and also for LT by Lemma 5.3 



It remains to verify the constraint corresponding to edges uv with u G S, v G V \ S. The 
contribution of Y^TeT-uveSfT) T\s^d ^(T) is unchanged. The following claim completes the proof of 
optimality. □ 

Claim 5.10. ^s( u ) + ^{S) < ^s( u ) + ^(S) with equality whenever uv G supp(x). 

Proof. 



(T) - tf(T) 



A(n(r)-*(r)) if res, 

-AA if T = S. 



Thus, 

V s (u) + *(5) = A(n s (u) - * s («)) + * (S) - + *s(«) + *(5) 
= X(U s (u) - - A) + + 

Now, ILj(it) — ^s(n) < A, and equality holds whenever uv G supp(x) n 5(S) by the consistency of 
and LT inside 5 (Lemma 5.1). □ 



Corollary 5.11. Let J 7 be a critical family with IT being an T-critical dual feasible solution. Al- 
gorithm Positively- critical- dual- opt in Figure^ transforms an arbitrary dual optimal solution ^ to 
an T- positively -critical dual optimal solution in at most \ J-\ iterations. 
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Proof. The correctness of the algorithm follows by Lemma 5.9 We bound the running time by 
showing that no set S G J- is processed twice. After a set S is processed, by Lemma 5.9, either 
IT and \P will be identical inside S or ^f(S) = 0. Once II and ^ become identical inside a set, it 
remains so during all later iterations. 

The value ^(S) could be changed later only if we process a set S' ~2 S after processing S. Let S' 
be the first such set. At the iteration when S was processed, by the maximal choice it follows that 
^(S') = 0. Hence fy(S') could become positive only if the algorithm had processed a set Z D S' , 
Z £ T between processing S and S", a contradiction to the choice of S' . □ 



5.3 Extremal dual solutions 

In this section, we prove Lemma [3. 2[ Assume J- C O is a critical family, with IT being an .F-critical 



dual. Let x be the unique optimal solution to Pj^(G,c) Let T x = {S € J- : x(S(S)) = 1} the 
collection of tight sets for x. A LT-extremal dual can be found by solving the following LP. 



min/i(^,n) 



E 



r(S) 
\S\ 



-r(S) < ^(S) - n(5) < r(S) VSgVUJ, 

^(S) = c(uv) Muv 6 supp(x) 

^(5) < c(in?) Vuu € E \ supp(x) 

> o vseT x 

The support of VP is restricted to sets in V D J- x . Primal-dual slackness implies that the feasible 



solutions to this program coincide with the optimal solutions of Djr(G, c) , hence an optimal solution 
to 



is also an optimal solution to Djr(G,c 



Lemma 5.12. Let T C O be a critical family with LT being an T -critical dual. Then, a U-extremal 
dual is also an J 7 -positively- critical dual optimal solution. 

Proof. We will show that whenever ^(S) > 0, the solutions ^ and LT are identical inside S. 
Assume for a contradiction that this is not true for some S £ J 7 . Let A = min{l, ^(5) /An,* (5*)} 
if An,* (5 ) > and A = 1 if A n *(S) = 0. Define ^ as in g. By Lemma 



DAG,c) 



and thus feasible to 



D 



5.9 



^ is also optimal to 



For every T € V U T x , let r(T) 

h(#,n) - h(v,n 



We show h(ty,Tl) < h(^,U), which is a contradiction. 
*(T) - n(T)| - |^(T) - n(T)|. With this notation, 



E 



r(T) 



If T \ S = 0, then ¥(T) = *(T) and thus r(T) = 0. If T C 5, T G _V U J 7 , then |tf (T) - U(T)\ = 
(1 - A)|^(T) - n(T)|, and thus r(T) = A|*(T) - n(T)|. Since tf(S) = tf(S) - AA, we have 
r(5) > -AA. 
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Let us fix an arbitrary u G S, and let 7 = max*rcs , :ueT,TeVu.F :E I^T 

r(T) 



h(^,U)-h(^,U)= 



2 E 

TQ&iteiyreVUF., 
A 



r(T) . r(5) 
|T| + |5| 



> 



7 



E 



|*(r)-n(r)| 



Tc5:«er,Tevu^ 



AA 

W\ 



> (n 5 ( u )-* s «) -_. 
7 p| 

Case i: If A > 0, then pick u £ S satisfying 115(7/) — \l/s(ii) = A. Then the above inequalities give 

1 1 



-/»(¥,!!) > AA( --p ) >U. 



The last inequality follows since IS"! > 7. 
Case If A = 0, then A = 1 and therefore, 

/»(*,n)-/i(*,n) > 1 



7 



E 



|¥(T) -II(T)| 



TcsiueT.Tevu^ 



Now, if II and \l/ are not identical inside S, then there exists a node u £ S for which the RHS is 
strictly positive. Thus, in both cases, we get h(^,H) < h(^f,JT), a contradiction to the optimality 
of * to \DU □ 



Proof of Lemma 3.2. By Lemma 3.1 the unique optimal x to Pj^(G,c) is proper-half-integral. 
Lemma 5.12 already shows that a T-extremal dual solution II is also ^-positively-critical. We 



need to show that the next family of cuts is a critical family. Recall that the set of cuts for the 
next round is defined as %' U H", where %' = {T G T : H(T) > 0}, and H" is defined based on 
some cycles in supp(x). We need to show that every set of H' U H" is II- factor-critical. This is 
straightforward for sets of T-C by the definition of the ^-positively-critical property. 

It remains to show that the sets of %" are also II-factor-critical. These are defined for odd 
cycles C G supp(x). Now, C G %" is the union of V(C) and the maximal sets Si,..., St of %' 
intersecting V(C). We have H(Sj) > for each j = 1, . . . ,£ and hence x(5(Sj)) = 1. 

Let u G C be an arbitrary node; we will construct the LT-critical matching M u . Let us contract 
all sets Si, . . . , Se to nodes si, . . . ,S£ w.r.t. II. We know by Lemma 4.1 iii) that the image x' of 
x is proper-half-integral and that the cycle C projects to a cycle in x' . The proper-half- integral 
property guarantees that C is contracted into an odd cycle C in supp(x'). Let v! be the image of 
u; the pre-image of every edge in C is a tight edge w.r.t. II in the original graph since II is an 
optimum to the dual problem Djr(G, c) Since C is an odd cycle, there is a perfect matching M', 



that covers every node in V(C) \ {u 1 } using the edges in C. 

Assume first u G Sj for some 1 < j < I. Then v! = Sj. The pre-image M of M', in the original 
graph contains exactly one edge entering each Sk for k 7^ j and no edges entering Sj ■ Consider the 
LT-critical matching M u for u in Sj. For k 7^ j, if a^b^ G M n 5(Sk), ctf. G S^, then, let M ak be the 
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II-critical matching for in S^. The union of M, M u and the M afe 's give a II-critical matching for 
u inside C. 

If u G C \ (uj =1 5j), then similarly there is a II-critical matching M afc inside every Sk- The 
union of M and the M afe 's give the II-critical matching for u inside C. We also have 11(5) > for 
all non-maximal sets S G W U "H" since the only sets with 11(5) = are those in T~L" , and they are 
all maximal ones. □ 



6 Convergence 



The goal of this section is to prove Lemmas 3.3 and 3.4 The first shows that the number of 



odd cycles in the support is nonincreasing, and the second shows that in a sequence of iterations 
where the number of cycles does not decrease, all the new cuts added continue to be included 
in subsequent iterations (till the number of cycles decreases). These structural properties are 
established as follows. First, we give an extension of Edmonds' primal-dual algorithm to half- 
integral matchings. Next, we argue that applying this algorithm to the current primal/dual solution 
leads to an optimal solution of the next LP. The analysis of the algorithm shows that the number 
of odd cycles is nonincreasing. Finally, we prove the extremal dual solution of the next LP must 
be the same as the one found by this combinatorial algorithm, and therefore the LP solution we 
use also satisfies the required properties. 

As mentioned in the introduction, our extension of Edmonds' algorithm to the half-integral 
setting and its application are for the purpose of analysis. We describe this algorithm first, derive 
its properties, then use it in the main proof of convergence of the cutting plane algorithm. 



6.1 The half-integral matching algorithm 

The next algorithm will be applied in certain contractions of G, but here we present it for a general 
graph G = (V, E) and cost c. We use the terminology of Edmonds' weighted matching algorithm 
[8] as described by Schrijver [231 Vol A, Chapter 26]. 

Let W C V, and let J- C O be a laminar family of odd sets that are disjoint from W . Let V w 
denote the set of one-element subsets of V\ W. The following primal 



(G,c) 



and dual 



Dj 



(G,c] 



arc 



Pf(G,c)) 



programs describe fractional matchings that leave the set of nodes in W exposed (unmatched) while 
satisfying the blossom inequalities corresponding to a laminar family T . The primal program is 
identical to Pj{G \ W, c) while optimal solutions to Djr(G \ W, c) that are feasible to Dj-(G, c) 
also optimal solutions to D^(G, c). 

min c{uv)x(uv) 

uv£E 

x(5(u)) = l VueV-W 
x(8(u)) = \/u G W 
x(8{S)) > 1 V5 G T 
x > 



(D${G,c)) 



max H-(S) 

S6V w UJ":ni)e5(S) 

n(5) > V5 G F 
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The algorithm is iterative. In each iteration, it maintains a set T C W, a subset C C T of cuts, 
a proper- half-integral optimal solution z to Pj(G, c), and an ^-critical dual optimal solution A to 
D^(G, c) such that A(5) > for every S G £. In the beginning T = W, C = J- and the algorithm 
terminates when T = 0. 

We work on the graph G* = (V*,i?*), obtained the following way from G: We first remove 
every edge in E that is not tight w.r.t. A, and then contract all maximal sets of C w.r.t. A. The 
node set of V* is identified with the pre- images. Let c* denote the contracted cost function and z* 
the image of z. Since E* consists only of tight edges, A(it) +A(v) = c*(uv) for every edge uv G E* . 
Since T is disjoint from W, the nodes in C will always have degree 1 in z*. 

In the course of the algorithm, we may decrease A(S) to for a maximal set S of C In this 
case, we remove S from C and modify G* , c* and z* accordingly. This operation will be referred 
as 'unshrinking' S. New sets will never be added to C. 

The algorithm works by modifying the solution z* and the dual solution A* . An edge uv G E* is 
called a O-edge/^-edge/ 1-edge according to the value z*(uv). A modification of z* in G* naturally 
extends to a modification of z in G. Indeed, if S G C* is a shrunk node, and z* is modified so 
that there is an 1-edge incident to S in G* , then let u\V\ be the pre-image of this edge in G, with 
u\ G S. Then modify z inside S to be identical with the A-critical-matching M Ul inside S. If there 
are two half-edges incident to S in G* , then let u\V\, U2V2 be the pre-image of these edges in G, with 
ui,U2 G S. Then modify z inside S to be identical with the convex combination (l/2)(M Ul + M U2 ) 
of the A-critical-matchings M Ul and M U2 inside S. 

A walk P = V0V1V2 ■ ■ ■ Vk in G* is called an alternating walk, if every odd edge is a 0-edge and 
every even edge is a 1-edge. If every node occurs in P at most once, it is called an alternating path. 
By alternating along the path P, we mean modifying z*(viVi+i) to 1 — z*(viVi+i) on every edge of 
P. If k is odd, vo = vp. and no other node occurs twice, then P is called a blossom with base vq. 
The following claim is straightforward. 

Claim 6.1 ([23, Thm 24.3]). Let P = vqVi . . .V2k+i be an alternating walk. Either P is an alter- 
nating path, or it contains a blossom C and an even alternating path from vq to the base of the 
blossom. □ 

The algorithm is described in Figure [7| The scenarios in Case I are illustrated in Figure |4} In 
Case II, we observe that T G B + and further, B + T\B~ = (otherwise, there exists a T—T alternating 
walk and hence we should be in case I). The correctness of the output follows immediately due to 
complementary slackness. We show the termination of the algorithm along very similar lines as the 
proof of termination of Edmonds' algorithm. 
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Half-integral Matching 

Input. A subset W C V, a critical family T C with all sets in .F disjoint from W, a 



proper-half-integral optimal solution w to (G, c) , and an F-critical dual optimal solution 



r to 

Output. A proper-half-integral optimal solution z to Pc(G, c) and an ^-critical dual optimal 
solution A to Pc(G,c) for some C C T. 

1. Initialize z = w, C = J 7 , A = V, and T = W. Let G* = (V*,E*), where E* C\ E are 
edges that are tight w.r.t. A, and all maximal sets of C w.r.t. A are contracted; c* 
and z* are defined by the contraction. Let R D T be the set of exposed nodes and 
nodes incident to ^-edges in z*. 

2. While T is not empty, 

Case I: There exists an alternating T-R-walk in G* . Let P = vq . . .V2k+i denote a 
shortest such walk. 

(a) If P is an alternating path, and v 2k+i £ T, then change z by alternating along P. 

(b) If P is an alternating path, and t>2fc+i £ R — T, then let C denote the odd cycle 
containing V2k+i- Change z by alternating along P, and replacing z on C by a 
blossom with base V2k+i- 



(c) If P is not a path, then by Claim 6.1 it contains an even alternating path P\ to a 
blossom C. Change z by alternating along P\, and setting z*(uv) = 1/2 on every 
edge of C. 

Case II: There exists no alternating T-R-walk in G* . Define 

B + := {S G V* : 3 an even alternating path from T to S}, 
B~ := {S £ V* : 3 an odd alternating path from T to S}. 



For some e > 0, reset 



A(S): 



(A(S) + £ ifS££ H 



\A(5)-e ifSefi-. 
Choose e to be the maximum value such that A remains feasible. 

(a) If some new edge becomes tight, then E* is extended. 

(b) If A(S) = for some S £ C n i3~ after the modification, then unshrink the node 
S. Set C:=C\S. 



Figure 7: The Half-integral Matching Algorithm 



Before we proceed to analyze the above algorithm, we briefly digress to note that it can be 
extended to a complete algorithm to find a minimum-cost perfect matching where the intermediate 
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solutions are half-integral and satisfy the degree constraints for all vertices. We next describe this 
combinatorial half-integral algorithm and then return to our main analysis. 

In the combinatorial half-integral algorithm, we make one modication to the Half-integral 
Matching subroutine. Namely, in Case 11(b) we keep track of the dual values on vertices in W 
and if one of them decreases to zero, then the corresponding vertex is uncontracted and the vertex 
is removed from W. 

The overall algorithm starts with the empty solution w = 0,T = 0,W = V, J- = Q and applies 
the following iteration repeatedly. 

1. [z,A,£] = Half-integral Matching(W, F, w,T). 

2. Contract cycles in z w.r.t. A; let the new W be the set of contracted vertices. 

3. For each vertex of the original W that is incident to a single edge in z, uncontract the vertex, 
extend z and A using the uncontracted edges and earlier dual values, and recurse this way 
till there are no contracted vertices incident to matching edges in z. Set the new primal and 
dual solutions to be the ones at the end of this procedure. Include each uncontracted set in 
the new T . 



We can prove the following about this algorithm. 

Theorem 6.2. Combinatorial Half-Integral is a strongly polynomial algorithm for minimim-cost 
perfect matching. 

The proof builds on that of the main subroutine and is deferred to the end of this section. 
Returning to the analysis of the Half-integral Matching algorithm, let f3(z) denote the number 
of exposed nodes plus the number of cycles in supp(z). We first note that /3(z) = /3(z*). This 



can be derived from Lemma 4.1 'iii) (We apply this Lemma in G \ T, observing that Pj(G,c) is 



identical to Pc(G\T,c)). Our next lemma (Lemma 6.3) shows that j3(z) is non-increasing. If (3{z) 
is unchanged during a certain number of iterations of the algorithm, we say that these iterations 
form a non- decreasing phase. We say that the algorithm itself is non-decreasing, if (3{z) does not 
decrease anytime. In the next section, we investigate properties of non-decreasing phases. These 
results will also show that every non-decreasing phase may contain at most | V| + \ J- \ iterations and 
therefore the algorithm terminates in strongly polynomial time. 

Lemma 6.3. Let z be an arbitrary solution during the algorithm, and let a be the number of odd 
cycles in supp(w) that are absent in supp(z). Then \ W\ + odd(w) > (3(z) + 2a. At termination, 
\W\ + odd(w) > odd(z) + 2a. 

Proof. Initially, /3(z) = \W\ +odd(w). Let us check the possible changes in /3(z) during an iteration 
of the algorithm. In Case 1(a), the number exposed nodes decreases by two. In Case 1(b), both the 
number of exposed nodes and the number of cycles decrease by one. In Case 1(c), the number of 
exposed nodes decreases by one, but we obtain a new odd cycle, hence /3(z) remains unchanged. 
In Case II, z is not modified. 

The only way to remove a cycle from supp(z) is by performing the operation in Case 1(b). This 
must be executed a times, therefore /3(z) < /3(w) — 2a. Further, there are no exposed nodes at the 
end of . Thus, on termination fi(z) = odd(z), and the claim follows. □ 
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6.1.1 The non-decreasing scenario 

Let us now analyze the first non-decreasing phase V of the algorithm, starting from the input w. 
These results will also be valid for later non-decreasing phases as well. In this case, we say that the 
algorithm is non-decreasing. Consider an intermediate iteration with z, A being the solutions, C 
being the laminar family and T being the exposed nodes. Recall that R 5 T is the set of exposed 
nodes and the node sets of the 1/2-cycles. Let us define the set of outer/inner nodes of G* as 
those having even/odd length alternating walk from R in G* . Let M and Mi denote their sets, 
respectively. Clearly, B + C M , B~~ C Mi in Case II of the algorithm. 

Lemma 6.4. If V is a non- decreasing phase, then if a node in V* is outer in any iteration of phase 
V, it remains a node in V* and an outer node in every later iteration of V . If a node is inner 
in any iteration of V , then in any later iteration of V , it is either an inner node, or it has been 
unshrunk in an intermediate iteration. 

Proof. Since V is a non-decreasing phase, Cases 1(a) and (b) can never be performed. We show 
that the claimed properties are maintained during an iteration. 

In Case 1(c), a new odd cycle C is created, and thus C is added to R. Let P\ =vq. . . V21 denote 
the even alternating path with vq £ T, V21 6 C. If a node u 6 V* had an even/odd alternating 
walk from vq before changing the solution, it will have an even/odd walk alternating from V2e £ R 
after changing the solution. 

In Case II, the alternating paths from T to the nodes in B~ and B + are maintained when the 
duals are changed. The only nontrivial case is when a set S is unshrunk; then all inner and outer 
nodes maintain their inner and outer property by the following: if u\V\ is a 1-edge and U2V2 is a 
0-edge entering S after unshrinking, with u\,U2 £ S, we claim that there exists an even alternating 
path inside S from u\ to U2 using only tight edges wrt A. Indeed, during the unshrinking, we 
modify z to M Ul inside S. Also, by the A-factor-critical property, all edges of M U2 are tight w.r.t. 
A. Hence the symmetric difference of M Ul and M U2 contains an alternating path from u\ to 112. 

We have to check that vertices in M — B + and Mi — B~ also maintain their outer and inner 
property. These are the nodes having even/odd alternating paths from an odd cycle, but not from 
exposed nodes. The nodes in these paths are disjoint from B~ U B + and are thus maintained. 
Indeed, if (B~ n M ) \ £> + / or (B + n Mi) \ B~ / 0, then we would get an alternating walk from 
T to an odd cycle, giving the forbidden Case 1(b). □ 

The termination of the algorithm is guaranteed by the following simple corollary. 

Corollary 6.5. The non- decreasing phase V may consist of at most \ V\ + (J 7 ! iterations. 

Proof. Case I may occur at most \W\ times as it decreases the number of exposed nodes. In Case 



II, either Mi is extended, or a set is unshrunk. By Lemma 6.4, the first scenario may occur at most 



\V\ times and the second at most times. □ 

In the rest of the section, we focus on the case when the entire algorithm is non-decreasing. 

Lemma 6.6. Assume the half-integral matching algorithm is non- decreasing. Let T be the initial 
dual and z, A be the terminating solution and C be the terminating laminar family. Let M and Mi 
denote the final sets of outer and inner nodes in G* . 

• If A(S) > T(S) then S is an outer node in V*. 
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If A(S) < T(S), then either S £ J- \ C, (that is, S was unshrunk during the algorithm and 
A(5) = 0) or S is an inner node in V* , or S is a node in V* incident to an odd cycle in 
supp(z) . 



Proof. If A(5) > r(5), then 5 £ B + in some iteration of the algorithm. By Lemma 6.4 this 



remains an outer node in all later iterations. The conclusion follows similarly for A(5) < T(S). □ 

Lemma 6.7. Assume the half-integral matching algorithm is non- decreasing. Let z, A be the 
terminating solution, C be the terminating laminar family and G* the corresponding contracted 
graph, M and M% be the sets of outer and inner nodes. Let : V* — > K be an arbitrary optimal 
solution to the dual Dq(G* , c*) of the bipartite relaxation. If S £ V* is incident to an odd cycle 
in supp(z), then A(S) = 0(5). Further S £ M implies A(5) < 0(5), and S £ Mi implies 
A(S) > 9(5). 

Proof. For 5 £ M UjVj, let £(S) be the length of the shortest alternating path. The proof is by 
induction on £{S). Recall that there are no exposed nodes in z, hence £(S) = means that 5 is 



contained in an odd cycle C. Then 0(5) = A(S) is a consequence of Lemma 5.8 both G and A 
are optimal dual solutions in G* , and an odd cycle in the support of the primal optimum z is both 
A-factor-critical and 0-factor-critical. 

For the induction step, assume the claim for £(S) < i. Consider a node U £ V* with £(U) = i + 1. 
There must be an edge / in E* between 5 and U for some 5 with £(S) = i. This is a 0-edge if i is 
even and a 1-edge if i is odd. 

Assume first i is even. By induction, A(S) < 0(5). The edge / is tight for A, and @(S)+Q(U) < 
c*(f). Consequently, A (U) > 8 (U) follows. Next, assume i is odd. ThenA(iS) > Q(S) by induction. 
Then, A(U) < Q(U) follows as / is tight for both A and 0. □ 

6.2 Proof of convergence 

Let us consider two consecutive solutions in Algorithm C-P-Matching. Let x be the unique proper- 
half-integral optimal solution to \Pj^(G, c)| and n be an ^-positively-critical dual optimal solution 
to Djr(G, c) We define ri! = {S : S £ T , 11(5) > 0} and %" based on odd cycles in x, and use the 
critical family % = H' U ri" for the next iteration. Let y be the unique proper-half-integral optimal 
solution to Pfi(G,c), and let ^ be an ^-positively-critical dual optimal solution to Du(G,c). We 



already know that II is an ^-critical dual feasible solution to D%(G,c) by Lemma 3.2 



Let us now contract all maximal sets 5 £ H with ^(5) > w.r.t. Vl/ to obtain the graph 



G = (V, E) with cost c. Note that by Lemma 5.8 II and are identical inside 5, hence this is the 



same as contracting w.r.t. II. Let x, y, II, and 'F be the images of x, y, II, and ^, respectively. 

Let H" = {S : 5 G H", (5) > 0}, and let W = WH" denote the union of the members of H". 
Let W denote the image of W. Then W is the set of exposed nodes for x in G, whereas the image of 
every set in U"\H" is an odd cycle in x. Let M = {T £ W : TnW = 0}, K = {T £ M : *(T) = 0} 
and M and JC be their respective images. All member of M \ K are contracted to single nodes in 
G; observe that K, is precisely the set of all sets in M of size at least 3. 

Using the notation above, we will couple the solutions of the half-integral matching algorithm 
and the solutions of our cutting plane algorithm. For this, we will start the Half-int egra l Matching 



algorithm in G with W, from the initial primal and dual solutions x and fl. Claim 6.8 n) justifies 
the validity of this input choice for . 
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Claim 6.8. (i) For every C C JC, y is the unique optimal solution to P^(G,c) and ^ is an 
optimal solution to D*(G,c). 

(ii) x is a proper-half-integral optimal solution to P^(G,c) and ft is a JC-positively-critical dual 
optimal solution to D^(G,c). 



Proof, (i) For C = 0, both conclusions follow by Corollary 4.3 - y is a unique optimal solution to 
Pu(G,c) and ^ is an %-positively-critical dual optimum. 

For an arbitrary C C JC, since y(6(S)) > 1 for every S G %' in the pre-image of C, y is a feasible 
solution to P^{G,c). Since y is optimum to the bipartite relaxation, this implies optimality of y 
to P^(G,c). Now, \l/ is optimum to D^{G,c) since satisfies complementary slackness with y. 
Uniqueness follows since a different optimal solution to P^{G, c) would also be optimal to Pq(G, c), 
but y is known to be the unique optimum to Pq(G, c). 



(ii) We will use Lemma 4.1 on an appropriate graph. We first setup the parameters that we 
will be using to apply this Lemma. First observe that x is the unique optimum to Pj^{G, c) as well 
as to Pu>(G, c) and II is an optimum to Du*{G, c). 

Let G w denote the graph G \ W, c w denote the cost function restricted to this graph. Let 
Tl w denote II restricted to the set M and x w denote the solution x restricted on G w . Since 
x(S(W,V \ W)) = 0, x w is the unique optimal solution to Px(G w ,c w ). By complementary 
slackness, H w is optimum to Dj\f(G , c ). 

We now apply Lemma |4.1[ iii) by considering the graph G w with cost c w , and the laminar odd 
family J\f. As noted above, x w and Vl^ are primal and dual optimal solutions. Let S G J\f\K, that 
is, y(S) > 0. Since S G M C H', we also have 11(5) = U W {S) > and hence S is a Il^-factor- 
critical set. Let us contract all such sets S w.r.t. U w . Then, by the conclusion of the Lemma, we 
have that x w is the unique proper-half-integral optimum to Pj^(G w , d w ). Since K is the set of 
nonsingular elements of N, it is equivalent to saying that x, is the unique optimum to P^{G,c). 

By Lemma 4.1 ii), we also have that f[ w is optimum to D^(G W , c w ) = D^(G \ W, c). Recall that 

Win t,\ \c „i„o +^ nlf/ u , 

W 



every optimal solution to D^(G \ W,c) that is feasible to D^(G,c) is also optimal to D^(G,c 
Observe that ft is on every node and subset of W. Consequently, all nonzero values of ft and ft 



^(G, c), we get that ft is also optimum to 
is /C-positively-critical by definition of JC. □ 



coincide. Since ft is feasible to D^(G : c), we get that ft is also optimum to D^{G, c). Further, ft 



Lemma 6.9. Suppose we start the Half-integral Primal-Dual algorithm in G, c, JC, W , from the 
initial primal and dual solutions x and II. Then the output z of the algorithm is equal to y. 



Proof. Since z is the output of the I 
P^(G,c) for some C C AC. By Claim 



alf-integral Matching algorithm, it is an optimal solution to 
i), y is the unique optimal solution to this program, and 
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consequently, z = y. □ 



Proof of Lemma 3.3. Using the notation above, let us start the Half-integral Matching algorithm 
in G, c, JC, W, from the initial primal and dual solutions x and II. Let z be the output of the 
half-integral matching algorithm. 



By Lemma 6.9, z = y. We first observe that odd(x) = \W\ + odd(x). This easily follows by 



Lemma 4.1 in), applied in G\W. Let a = \H"\H"\. There is an odd cycle in supp(rr) corresponding 
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to each set of %" \ %" . None of these cycles may be contained in supp(z) = supp(y) as otherwise 
the corresponding cut in H" would be violated by y. Thus Lemma |6.3| implies odd(y) = odd(z) < 
\W\ +odd(x) — 2a and Lemma [47T](iii) implies odd(y) = odd(y). Hence, odd(y) < odd(x) — 2a. □ 



The next claim formulates a simple consequence of the above proof. 

Claim 6.10. If odd(y) = odd{x), then H" = %" . Further, the Half -integral Matching algorithm 
applied in G, c, tC, W , with starting solution x, LI is non- decreasing. □ 



This claim already implies Lemma 3.4 for j = i + 1. To analyze the scenario odd(x) = odd(y), 
the results in Section 6.1.1 are applicable since the algorithm was non-decreasing. Let us start the 
half-integral matching algorithm in G, c, JC, W, from the initial primal and dual solutions x and 
fi. Consider the final dual solution A with corresponding laminar family £ and define A in G as 
follows. 

If S C T for some T G %, Vl/(T) > 0, then set A(S) = ^(S) (this defines the dual solutions for 
sets and nodes inside T that were contracted to obtain G). If S G £ U V, then set A(S) = A(S) 
for its pre-image S (this defines the dual solutions for sets and nodes on or outside T that were 
contracted to obtain G). 

Claim 6.11. A is an optimal solution to D-}{(G,c). 

Proof. We show feasibility here. Optimality follows by complementary slackness with y and the 
definition of shrinking. The nontrivial part to show feasibility is to verify that A (5) > if S E H. 
This is straightforward if S Q T for some T & 7i, ^f(T) > 0, or if S G £■ Let S G % be a maximal 
set with ^(S) > 0, and let s denote its image in G. Then during the Half-integral Matching 
algorithm, A (s) m ight have decreased below 0. We show that this cannot be the case. 



By Claim 6.8 1), \l/ is an optimal solution to the bipartite relaxation Dq(G, c). If A(s) < 0, then 
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A(s) < < n(s). We started the algorithm with the dual solution L = n, therefore by Lemma 
either (1) s is an inn er node or (2) s G V* and s is incident to an odd cycle in supp(z). In both 



cases, by Lemma 6.7, we get that A(5) = A(s) > = *$>(S) > 0. □ 



Lemma 6.12. Assume odd(x) = odd(y) for the consecutive solutions x and y. Then A = ^ and 
hence A = 



Proof. By Lemma 5.8 if S is such that ^(S) > 0, then ^ is identical to II inside S, and so by 



definition, also to A. Therefore, it suffices to prove A = ^. By Claim 6.8 i), \E' is an optimal 
solution to Dq(G, c). The following claim completes the proof by showing that h(A, H) < h(^f, II), 
and equality can hold only if A and \& are identical. □ 

Claim 6.13. For every S G VU/C, |A(5) - fl(S)\ < \$(S) - fl{S)\ and equality holds only if 
A(5) = 

Proof. The claim will follow by sh owing that for every S G V U fC, either H(S) < A(S) < ^(S) 



or n(S) > A(S) > ^(S). By Claim 6.10, odd(x) = odd(y) implies that the half-integral matching 
algorithm was non-decreasing. Let z be the terminating primal solution, and let G* = (V*,E*), c* 
be the corresponding contraction of maximal sets of £ in G, c. Note that ^f(S) = if S G £. Let 
G : V* ->■ M be defined as follows: 



9(5* 



tt(S)-A A> ft(S) ifSGA 
$(S) if S G V* \ £. 
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Claim 6.14. is an optimal solution to Dq(G* ,c*). Further, @(S) < holds for every Sg£. 
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5.1 



, y and ^ are optimal 
is consistent with A 



Proof. If S G £, then A(5) > and hence y(5(S)) = 1. £ is a pro per odd family in G with A 
being an ^-critical dual (by definition of the algorithm). By Claim 
solutions to Px(G,c) and to D^(G,c), respectively. Hence, by Lemma 
inside every set S G C such that y(S(S)) = 1. Further, A^^(S) > for every S G C. Thus, by 

Lemma 5.2 is an optimum dual solution to Dq{G*,c*). The second part follows since ^(S) = 
and A A ^(5) > if S G C. □ 



Now consider S C V. If A (S) > 11(5), then by Lemma 6.6, we have that S G V* and S G M Q 



Consequently, by Lemma 6.7, ®{S) > A(S). If S G £, then > ©(5*) > A(5), a contradiction to 



A(5) > 0. Thus, S = {s} for some s G V. Hence, *(s) = 9(s) > A(s) > n(s 



If A(5) < n(5), then by Lemma 
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we have that either (1) S G /C \ C, that is, A(5) = and 
S was unshrunk or (2) S G A/i or (3) S G V* and 5 is incident to an o dd cycle C in supp(z). If fij, 
then #(5) = = k(S) < 11(5). In both cases (%) and ^j, Lemma [6/7 
then ^(5) = < A(5) < If 5 = {s} for some s G V", then *(s) 



gives 9(5) < A(5). If 5 G £, 
0(a) < A(s) < n(s). □ 

Proof of Lemma \3.4\ Let Xi be the solution in the i'th iteration (above, we used x = Xi and 
y = Xi + \). Assume the number of odd cycles does not decrease between iterations i and j. By 
Claim [6?10 if we run the half-integral matching algorithm between Xk and Xk+i, for i < k < j, it 



is always non-decreasing. 

We first run on the co ntra cted grap h G 
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and 



solution n = Hj. Lemmas 
y = Xi + i and dual optimal solution A 
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Gi starting from primal solution x = Xi and dual 



For j = i + 1, the statement follows by Claim 



show that it terminates with the primal optimal solution 


since %" = %" means that all cuts added 
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in iteration i have positive dual value in iteration i + 1. Further, all sets in %" were contracted to 



exposed nodes in Xj. By Lemma 6.6, these will be outer nodes on termination of the half-integral 
matching algorithm as well. Let G* be the contracted graph upon termination of the Half-Integral 
Primal-Dual algorithm. 

Let J = J'UJ" be the set of cuts imposed in the (i+2)'th round, with J = {S G U : V(Z) > 0}, 
and let J" be defined according to odd cycles in xi + \. Let $ be the extremal dual optimal solution 
to Dj(G,c). 

Let us run the half-integral matching algorithm from Xj-fi to £j+2- We start the algorithm with 
the contracted graph Gj+i, which results by contracting all sets with $(5) > 0, S G J . Let G* +1 
be the initial contraction of Gj+i used by the algorithm. 

The key observation is that while the underlying graphs Gi and Gi+i are different, G* +1 can 
be obtained from G* by contracting those odd cycles corresponding to the sets of J" . Every other 
node that was inner or outer node in G* will also be inner or outer node in G* +1 , including the 



members of T~L" . By Lemma 6.6, the members oiH" will be outer nodes at termination, along with 
the new outer nodes J" . 

Iterating this argument one can show that every set that was imposed based on an odd cycle 
between iterations % and k will be outer nodes at the termination of the Half-integral Matching 
algorithm from x^ to Xfc+i. □ 

We conclude this section with a proof of convergence of the combinatorial half- integral algorithm 



(Theorem 6.2 ). 
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Proof of Theorem 6.2 We need to verify two things. First that the inputs to the main subroutine in 
each iteration are valid, and second that the algorithm makes progress towards an integral solution. 

To see the first, we simply use the contraction properties, namely any set that is uncontr acted 
has a half-integral cycle, which can then be adjusted to maintain a half-integral solution after 
uncontraction. Moreover, the dual values for vertices inside will be the values at the point of 
contraction. An argument similar to the proof of Lemma 6.9 establishes the validity of this dual 
solution. 

To show progress, we use odd(x) as in the analysis of the cutting plane algorithm. In the 
cutting plane analysis, in a non-decreasing phase, the cuts that are added in each iteration remain 
in all iterations until the end of the non-decreasing phase. Instead, what we argue here is that the 
number of contracted vertices plus the size of the laminar family at the end of a phase is at least 
the total number of half-integral cycles encountered during the course of the non-decreasing phase. 
To see this, we first note that each half-integral cycle at the end of a phase is contracted. Next, 
when a vertex is uncontracted, the corresponding set is added to the laminar family. Now we need 
to verify that when a set is added to the laminar family it remains in the family till the end of the 
phase. This is similar to the proof of Lemma 3.4 Therefore the length of a non-decreasing phase 
is bounded, and polynomiality of the algorithm follows. □ 



7 Uniqueness 



In Section |5.2| we proved that if J 7 is a critical family, then there always exists an J^-positively- 
critical optimal solution (Corollary |5.11[ ). This argument did not use uniqueness. Indeed, it will 
also be used to derive Lemma [2. 3[ showing that a perturbation of the original integer cost function 
satisfies our uniqueness assumption Q. We will need the following simple claim. 

Claim 7.1. For a graph G = (V, E), let a,b : E —> be two vectors on the edges with a(5(v)) = 
b(5(v)) for every v E V. If a and b are not identical, then there exists an even length closed walk 
C such that for every odd edge e £ C, a(e) > and for every even edge e £ C, 6(e) > 0. 

Proof. Due to the degree constraints, z = a — b satisfies z(6(y)) = for every v E V, and since a 
and b are not identical, z has nonzero components. If there is an edge uv E E with z(uv) > then 
there must be another edge uw E E with z(uw) < 0. This implies the existence of an alternating 
even closed walk C where for every odd edge e E C, < z(e) = a(e) — 6(e), and for every even 
edge, > z[e) = a(e) — 6(e). This proves the claim. □ 



Proof of Lemma \2.3[ Let c denote the perturbation of the integer cost c : E — > Z. Consider a 
graph G = (V,E), perturbed cost c and critical family T where Q does not hold. Choose a 
counterexample with \T\ minimal. Let x and y be two different optimal solutions to Pj^(G,c). 
Since J 7 is a critical family, by Corollary |5.11 there exists an ^-positively-critical dual optimal 



solution, say II. 

First, assume T = 0. Then x and y are both optimal solutions to the bipartite relaxation 



Pq(G,c). As they are not identical, Claim 7.1 gives an even closed walk C such that x(e) > on 
every even edge and y(e) > on every odd edge. Let 71 and 72 be the sum of edge costs on even 
and on odd edges of C, respectively. Then for some e > 0, we could modify x by decreasing x(e) 
by e on even edges and increasing on odd edges, and y the other way around. These give two other 
optimal matchings x and y, with c T x = c T x + (72 — 71 )e and c T y = c T y + (71 — 72)^- Since x and 
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y are both optimal, this gives 71 = 72- However, the fractional parts of 71 and 72 must be different 
according to the definition of the perturbation, giving a contradiction. 

If T / 0, after modifying x, y to non-identical optimal solutions a, b respectively (if necessary), 
we will similarly identify an alternating closed walk C in supp(a) U supp(6) with the additional 
property that if C intersects 5{S) for some set S G P, then it does so in exactly one even and one 
odd edge. The modifications a and b defined as above would again be feasible, implying 71 = 72. 

We first claim that H(S) > must hold for all S £ T . Indeed, if H(S) = for some S 6 P, then 
x and y would be two different optimal solutions to Pjr^ s ^(G, c), contradicting the minimal choice 
of P. Let T £ P U {V} be a smallest set with the property that x(5(u, V — T)) = y(5(u, V — T)) 
for every u £ T, but x and y are not identical inside T. Note that V trivially satisfies this property 
and hence such a set exists. Let S denote the collection of maximal sets S in P such that SCT 
(S could possibly be empty). 

Inside each maximal set S £ S, we modify x and y such that they are still both optimal and 
different from each other after the modification. Since H(S) > 0, we have x(5(S)) = 1, y(5(S)) = 1 
and S is II-factor-critical. For u £ S, let M u denote the II-critical matching for u inside S. Let 
a u := x(5(u,V - S)),a' u = y{5(u,V - S)) and w := J2ueS a ^ M u,w' := Y^ u es a 'u M u- If x and y 
are not identical inside S and a u = a' u for every u £ S, then S contradicts the minimality of T. 
Hence, for each maximal set S £ S we have two cases - either (1) x and y are identical inside S, 
and thus for every u £ S , a u = a' u or (2) a u 7^ a' u for some u £ S. 

Fix S £ S. The following modified solutions 




x(e) if e£5{S)UE[V\S], 
w{e) if e € E[S\, 



b(e) 



y(e) if e £ S(S)U E[V\S], 
w'{e) if e £ E[S], 



are also optimal solutions since we used only tight edges w.r.t. n for these modifications. 

Claim 7.2. Consider two edges u\Vi £ 5(S) n supp(a), U2V2 £ 5(S) n supp(b), u±,U2 £ S. Then 
there exists an even alternating path P$ inside S between u\ and U2 such that a(e) > for every 
even edge and 6(e) > for every odd edge. Also, consider the extended path P' s = u\V\PsU2V2- If 
there exists a set Z C S, Z £ T such that V(P' S ) D Z 7^ 0, then P' s intersects 5(Z) in exactly one 
even and one odd edge. 

Proof. By the modification, supp(a)ni?[S'] contains the n-critical-matching M Ul and supp(6)ni?[S'] 
contains the n-critical-matching M U2 . Then the symmetric difference of M Ul and M u , 2 contains an 



U\ — U2 alternating path satisfying the requirements (by Lemma 2.2). □ 



We perform the above modifications inside every S £ S, and denote by a and b the result of all 
these modifications. Let us now contract all sets in S w.r.t. n; let G',c',F denote the resulting 



graph, costs and laminar family respectively. By Lemma 4.1 ii) , the images a' and b' are both 
optimal solutions in P^,(G' ,c'). 

We claim that a' and b' are still not identical inside T' . Indeed, assume a' and b' are identical 
inside T' . Then we must have had case (1) for every S £ S. Note that we only modified x and y 
on the edge set U^gs-E^iS 1 ], and hence the contracted images must be identical: x' = a' and y' = 6', 
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and consequently, x' = y' . As x and y are not identical inside T, they must differ inside at least 
one S £ S. This contradicts the fact that case (1) applied for every S G S. 

By the assumption on T, we also have a'(5(u', V — T')) = b'(5(u' , V — T')) for every u' £ T' . 



Hence Claim 7.1 is applicable, giving an alternating closed walk C with a'{e) > on every even 
edge and b'(e) > on every odd edge. Now, we can extend C to an even alternating closed walk 
C in the original graph using the paths Ps as in the above claim. The resulting closed walk C will 
have the property that if there exists a set Z C S, Z E T such that V(C) f]Z ^ 0, then C intersects 
6{Z) in exactly one even and one odd edge. □ 

8 Open Questions 

Our initial motivation was to bound the number of iterations of the cutting plane method using 
the Padberg-Rao procedure. This question remains open and any analysis would have to deal with 
non-half-integral solutions. 



Within our algorithm, Lemma 3.1 shows that it is sufficient to use positively-critical dual 
optimal solutions to maintain proper-half-integrality. Can we prove efficient convergence of our 
cutting plane algorithm using positively-critical dual optimal solutions (without using extremal dual 
solutions)? We believe that such a proof of convergence should depend on whether the following 
adversarial variation of Edmonds' algorithm for perfect matching is polynomial time. Suppose 
we run the Edmonds' perfect matching algorithm, but after every few iterations, the adversary 
replaces the current dual solution with a different one, still satisfying complementary slackness 
with the (unchanged) primal solution. 

Given the encouraging results of this paper, it would be interesting to prove efficient convergence 
of the cutting plane method for other combinatorial polytopes. For example, one could try a similar 
approach for finding an optimal solution for 6-matchings. Another direction could be to try this 
approach for optimizing over the subtour elimination polytope. 
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